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abstract: 

If  an  elastic  bar  with  some  ihitiall  imperfection  is.  subjected  to 
increasing  compressive  axial!  load.,,  the-  lateral  deflection,  .measured 
anywhere  along  Its  length,,  in  creases;  monofonically;  the  primary  buck- 
ling mode  becoming  more:  and  more;  dbmihanlr:  in  the  deflected  shape..  .  A 
classical  approximate  technique:,,  due:  to  R\.  V. .  Southwell,  .correlates  : 
axial  load  and  lateral  deflection  for  determination  of: the  primary 
buckling  load  of  the.  bar..  However,.  Southwell' s  approximation  may  be 
inaccurate  if  the  primary  mode  does:  not:  predominate  in:  the  :def  lected 
shape.   The  present  thesis  proposes  a.  technique  which  employs  axial 
loading,  but  also  transverse  loading  to  insure  dominance  of  the 
primary  mode  in  the  deflected  shape. 

The  technique  allows  experimental  determination  of  the  second 
buckling  load,  again  using  appropriate  transverse  loads  to  insure 
dominance  of  that  mode  in  the  deflected  shape. 
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nomenclature: 

a  Coefficients  which  characterize  deflection- due  to  loads 
n 

b  Coefficients  which  characterize  initial  .deflection 
n 

e.  Error  associated  with  the  ith  data  point - 

1  —      r 

h  Height  of  column  cross-section  (in  y  direction) 

i,j,n  Integer/  subscripts: 

k  Lateral  stiffness  of  column  (lateral  load  to  produce 
unit  lateral  deflection) 

v  Deflection  divided  hy  corresponding  load  (6/P) 

w  Width  of  column  cross— section 

x  Axis  parallel  to  applied  axial  load 

y  Axis- parallel  to.  applied  lateral  load 

B  Flexural  rigidity  EI 

C  Compact  notation  for  specified  function  of  T  and  P 

n  ^                 r                         n 

E  Young's  modulus 

I  ,J  ,K  Constants  arising  from  evaluation  of  certain  specified 

n  n  n      .      , 
integrals 

L  Effective  length  of  column 

P  Compressive  axial  load 

P  Buckling  load  for  the  nth  mode 

n  — 

Q  Lateral  load 

S  Sum  of  the  squares  of  e. 

T  Tensile  axial  load 

U  Strain  energy 

V  Voltage 
W  Work 

Y  Orthogonal  functions 
n  to 


6  Measured  deflecrtxan 

p  1/k 
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I.   INTRODUCTION' 

The  problem  of  determining  the  elastic  buckling-  load  off  a-  column 

has  been  one  of  both  theoretical  interest:  and  practical  importance . . 

Euler's  theory  for  initially  straight  columns  has  been,  accepted:  and 

widely  used  for  many  years.   R-  V".  Southwell  Q5J  treated  the:  problem 

of  initially  imperfect  columns  from  both  theoretical  and:  experimental- 

viewpoints.   If,  before  loading,  the  column  has  some,  small"  distortion, 

usually  a  superposition  of  many  modes,  the  application  of '  a'-  compressive 

axial  load  P  magnifies  the  amplitude  of  the  nth  mode  by  a  factor. 

P  /(P  -P),  where  P   is  the  nth  mode  buckling  load,.  As:  P  increases, . 
n   n  n         & 

magnification  of  the  higher  mode  Imperfections  (n  >  1)  becomes  rela- 
tively insignificant  as  compared  to  that  of  the  first  mode.   If"  the 
column  is  prismatic  and  the  deflection  is  measured  at  mid-span,  even 
numbered  modes  make  no  contribution  at  all,  and  the  dominance  of  the 
first  mode  is  even  more  pronounced.   Southwell  postulated  that  to  a 
close  approximation  5  =  5  /(1-P/P  ),  where  <5   is  the  first  mode  imper- 
fection and  <5  is  the  mid-span  deflection.   Alternatively  P.v  -6  =  6 
where  v  =  6 /p.   Accordingly,  If  a  plot  of  v  vs.  <5  is  made,  a  straight 
line  with  inverse  slope  equal  to  P   results,  provided  P  is  large 
enough  to  assure  first  mode  dominance. 

Successful  use  of  Southwell's  method  for  determining  buckling 
loads  depends  upon  the  unloaded  column  having  an  initial  lateral 
deflection  in  which  there  is  a  significant  component  of  the  first 
buckling  mode.   If  this  is  not  the  case,  the  desired  linearity  of 
(an  appropriate)  plot  will  not  be  observed  and  it  will  not  be 
possible  to  make  an  accurate  determination  of  the  first  buckling 


load.   This  paper  presents,  an.  extension  off  Southwell's  methods.  The- 
proposed  method  does  not:  rely  cm.  accidental!  imperfection. .  Lateral- 
loads  are  used  to  insure  dominance  of :  the:-  nth  mode  in  the  deflected 
shape  when  determining  the  nth  buckling-  load..  By  using  appropriate 
lateral  loads,  the  buckling  load'  for  any  mode  may  be  determined,. 
although  only  the  first  two  are  considered: here . 
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II.      THEORETICAL   ANALYSIS 


■r 


y+ 


Eigure  2.1 
Loading  Diagram 


A.   GENERAL  CASE 

Eigure  2.1  represents  a  bar,  with  either  variable  or  constant 
cross- sect ion,  having  some  initial  imperfection  and  subjected  to 
axial  loading  P  and  lateral  loading  Q.   All  forces  and  axes  are 
shown  in  their  positive  sense.   The  solution  for  the  deflected  shape 
may  be  found  using  a  linear  combination  of  orthogonal  functions  to 
represent  the  initial  deflection  and  the  deflection  under  load.   The 
case  represented  here  is  for  a  single  transverse  load,  but  the  method 
may  be  extended  for  any  number  of  loads.   Here  and  elsewhere  it  is 
assumed  that  all  loads  and  deflections  are  confined  to  a  single  plane, 
namely  the  x-y  plane. 

Let  y  be  the  initial  ordinate,  measured  from  the  chord  AB  and  let 
y.  be  the  deflection  due  to  external  loads.   The  total  ordinate  after 
the  loads  have  been  applied  is  y  =  y  +  y  . 


Assume  y   =  b  Y  +  b  Y.  +  b.Y  +  ... 
o    11    22     33 


y±   =   alYl  +  a2Y2  +  a3Y3  + 


(2.1) 
(2.2) 
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where  Y      =  Y   (x)    (n  =  1,2,3,.......)    £s  a  sret:  off  orthogonal  functions: 

n  n 

in  the   interval   (0,L)  which  satisfy  the,  differential-  equation  B-lai 
subject   to   conditions  B-lh    (see  Appendix  B)..     As- the   loads   are 
applied,   the   distance  between  A  and  B   decreases   by  an  amount 


X   -   X 


L  2 

i  r    r  ^i^q*" 


dx 


dx  - 


Z 


n^y. 


dx 


dx 


-\\     [Cy')2  -CyD2]  dx 


(2:3> 


where  y'    =   ^(a.+b.)   YT;   y''   =  ^a.YT.;:  yr  =  ^b.Y.l'. 

J  ii       x       1  x.  r      a  ri 

and  the  prime  denotes   differentiation,  with,  respect,  to  x. 


(yD2    =   EZ(a.+b.)(a.+b  JY!'   YV.;:  (yD2=  £Xb  ..  h  ..  Y.V  Y'. 
i      i        ]      :      i      1  o  i    I    l    ] 

E 


;.    x  -  x    =  £  e 

o        2 


(a.+b.)2-  b.2 

.     ii  i  • 


(YV)zdx 


a.' 


=   Z    (~     +   a.b.)  (Yl)    dx 

2  ii/  l 


=  Z  (-i-     +   a.b.)    J 
2  l   l        n 


(2.4) 


For  definitions  and  properties  of  the  functions  appearing  here  and 

below,  see  Appendix  B. 

Now  consider  the  energy  variations  due  to  some  small  additional 

deflection  da  Y   and  assume  da   is  positive  in  the  same  sense  as  Q. 
n  n  n    r 

The  work  done  by  the  force  P  is 


3(X-  X  ) 

dW_  =  P  — 2_  -da   =  -  P(a  +b  )•  J  •  da 

P       oa         n      n  n    n    n 
n 


(2.5) 


The  work  done  by  the  force  Q  (located  at  x  =  x  )  is 

dWn  =  QY  (xn)-  da 
Q     n  0     n 


(2.6) 
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The  increase  in  strain  energy-  is 


dU  = 


3U_ 
3a 


da  = 


n 


n   3a 


n 


^TcTCx)  (y")2^ 


da 


n 


n   L  jj  J- 


where  B  =  El(x)  ::  B(x) 
3 


dU  = 


3a 


n 


J-  ir  F  a.,  a..  /  BE  YV.  Y\r  dxl  da 
2     ^  x  j  J     x  j   J-   rr. 


3   1     2 

3a   2    x  x 
n 


=  a  K  da 
n  n   n 


Q 
da 


n 


(2.7) 


Equating  the  increase  in  strain,  energy  and  the  work  done  by  the 
external  forces ,  an  equation  for  the  coefficients  a  may  be  found 


P(a  +b  )J  •  da  +  QY  (xn)  -da  =  a  K  da 


n  n  n    n    '  n   0' 

QY  (xn)  +  Pb  J 

'  n  0 n  n 

an  "    K  -  PJ 
n     n 


or  since  K  /J 

n  n 


a  = 
n 


§-  Y  (xn)  +  ^b 
K   n  u    P   n 
n  n 


^f- 


n 


y  =  v±  +  y0  =  £ 


K  Y   (xn)  +  b 

n  n   u n 

i-l- 


n  n   n 


Y  (x) 

n 


(2.8) 


(2.9) 


These  terms  may  be  ordered  so  that  P  <  PQ  <  P_  • • • 

_L     A.  o 

It  is  obvious,  for  values  of  P  approaching  P  ,  that  if  Q  is  sufficiently 
large,  the  first  term  of  the  series  is  strongly  dominant.   Southwell, 
rather  than  applying  a  lateral  load  Q,  depended  on  some  accidental 
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imperfection  such  that  h     £   Q-  Siaweswear,,  if:  th &.  nature  of  the:  initial- 
imperfection  is  such  that  h   >>  b  fear  same:  rr  »  11,.  the  firsts  term-  does  : 
not  dominate  unless  P  is  very  nearly- equal,  txr  F..  Consequently,  .the -:■ 
linear  relation  which  Southwell's  analysis  depends  on  may  not"  appear: 
in  plots  of  experimental  data. 

Equation  (2.9)  also  demonstralres  something:  which  does  not: seenr  to 
have  been  treated  in  previous  literature:..  If"  the"  value  of  ."the :  applied  : 
lateral  load  Q  is  sufficient,,  the  first  terra  will  dominate,  even,  if  :"bT  . 
is  not  the  largest  b   and  even,  if  P  daes~  not:  approach  P"  .   Furthermore, 
if  the  lateral  load  is  applied  near  the  column  midlength ,  then 
Y  (x  )/K   >>  Y  (x  )/K  .   Thus,,  in.  effect,,  the:  first  mode  imperfection 
may  be  increased  by  the  use  of  lateral  loading,,  insuring  dominance  of: 
the  first  mode  in  the  deflected  shape,,  regardless  of  the  value  of. 
P(P  <  P1).   By  increasing  the  effective  first  mode  imperfection  to 
b  *  =  b  +  (QY  (x  )/K  ),  the  accuracy  and  reliability  of  Southwell's 
method  can  be  greatly  enhanced. 

Another  experimental  procedure,  which  involves  observing  deflection 
as  a  function  of  lateral  load  Q  for  a  few  particular  values  of  P,  will 
be  described  subsequently. 

For  the  special  case  of  columns  which  are  symmetric  with  respect  to 

a  normal  plane  through  the  center,,  if  the  lateral  load  is  applied  at 

mid-span,  non-zero  values  of  the  coefficients  a  occur  only  for 

n  J 

n  =  1,3,5,  ...   These  coefficients  attenuate  rapidly  (l>>a_/a  >>ar/a,  .. 

o  1        o      1 

and  it  is  evident  that  the  first  term  of  equation  (2.9)  dominates  the 
others  even  if  1-P/P   is  not  small. 

This  suggests  that  the  value  of  P  might  be  accurately  approximated 
even  if  tensile  loadings  (i.e.,  negative  values  of  P)  are  used.   First 
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term  dominance  is  diminished  slightly,  but  not  sign ffrcanrly-, .  for 
tensile  loading.   In  the  experimental  program  to.  be  described,,  tensile  - 
loading  was  used  partly  for  convenience,  but  largely  to:  demonstrate 
what  appears  not  to  have  been  done  before,  namely  to  evaluate  a 
compressive  critical  load  without  ever  applying  compression.. 

Assuring  that  the  applied  axial  compressive  force  F  does  not  closely 
approach  the  first  buckling  load  P  permits  the  experimeTrtaX  determi- 
nation of  the  second  critical  buckling  load!!  To  do  this  it  is  necessary 
to  assure  dominance  of  the  second  term  of  equation  (Z..9)..  This  can  be 
done,  again  using  lateral  loads  to  "excite"  mainly  the  second  mode.. 
Generally,  it  is  difficult  to  develop  an.  experimental  technique:  which 
will  permit  accurate  establishment  of  the  true  first-  made  shape:  Y  (x) 
and  also  one  which  is  orthogonal  to  It.   However,,  if  the  column  is 
symmetric  with  respect  to  a  central  normal  plane,  clearly  the  mode 
shapes  are  alternately  symmetric  and  antisymmetric.   Thus,  antisym- 
metrical  loading  can  "excite"  only  the  second  and  higher  antisym- 
metrical  modes.   If  the  lateral  loading  consists  of  two  oppositely 
directed  forces  near  the  quarter  points,  fourth  mode  excitation  will 
be  small  and  the  second  should  dominate. 

B.   PRISMATIC  BAR 

Although  there  is  no  necessity  for  determining  the  buckling  load 

of  an  elastic,  prismatic  bar  by  experimental  means  since  the  theory 

provides  complete  information  for  that  case,  it  is  still  reasonable 

to  select  such  a  bar  for  experimental  confirmation  of  the  general 

analysis.   The  particular  form  of  the  preceding  results  for  the  case 

of  a  prismatic  bar  may  -be  easily  found.   The  orthogonal  functions  are 

Y   =  sin  n7Tx/L.   The  constants  K   and  J   are  then 
n  n      n 
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EC     = 

EL 


EInV 

a 

2L 


(2:io) 


2    2 
Jn   "  ~2L~ 


(2.11) 


Equation   (2.9)    is   then 


r  ,    2QL3  .      n7re0 

/  — -  sin 


y  =  y-l  +  y0  =  z 


ElnV 


+  b 
L      n 


sin. 


mrx 


V 


n 


j. 


(2:12) 


For  first  mode  determination,,  the  load.  Q  could.'  he:  applied  anywhere: in 

(0,L)  and  the  deflection  could  be  measured  at  the  same  or  any  other 

position  in  (0,L).   However,,  it  is  convenient  to.  do  both  at  the.  center,  . 

x  =  L/2.   If  this  is  done,  then  x  -   x  —  11/21  and  equation  (2.12)  becomes 

a 


y  = 


2QL 


EItt1 


x-\    81(1-9?7)    625a-hf1) 


+  . 


(2.13) 


1- 


9P. 


If  instead  of  a  compressive  load  P,  there  is  a  tensile  load  T,  clearly 
T  =  -P  and  equation  (2.13)  becomes 


y  =£ 

y       k 


1        1 
+  


T  T 

1  +  f-  81(1  +  ±-  ) 

1  1 


+•••  +  y 


(2.14) 


where  the  quantity  k  represents  lateral  stiffness 


k  =  EIttV2L3 


and 


y* 


1 


3    + 


T 


1  + 


9P. 


1  + 


25P. 
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With  the  application  of  Q  at  mid-span,  the  bar  will  deflect  essentially 
in  a  first  mode  configuration,  higher  mode  contributions  being  small. 

C.   ESTABLISHMENT  OF  REFERENCE  LEVEL 

The  last  term  in  equation  (2.14)  presents  a  problem  in  that  its 
exact  value  is  unknown  and  cannot  be  determined  in  the  laboratory.   The 
problem  may  be  circumvented  however.   Referring  to  Fig.  (2.2),  three 
load  cases  are  shown: 
(a)   T  =  Q  =  0; 


T  T 


Figure  (2.2) 
Experimental  Loading  Cases 

(b)   Q  =  0,  T  =  T;  and  (c)  Q  =  Q,  T  =  T.   With  no  loads  applied  there 
is  some  initial  imperfection  b.   When  load  T  is  applied  (Q  =  0),  the 
amplitude  of  the  imperfection  is  reduced  to  an  amount  y*  which  is 
exactly  equal  to  the  last  term  in  equation  (2.14).   With  the  load  T 
applied  and  maintained,  a  reference  level  may  be  established  for  all 
subsequent  deflections  due  to  lateral  loading.   Now  as  the  lateral 
load  Q  is  applied  (holding  T  constant),  the  deflection  measured  in 
the  laboratory  is  simply 

"ftf-V+    1  T    + — S— +---'         (2-15) 

\l+¥-       81(1  +T!F)    625(1  t^) 
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D.   NON-PRISMATIC  BARS 

For  non -prismatic  bars,,  equation  (2.3)   describes  the  deflected 
shape  In  terms  of  the  applied  loads,,  the  initial!  imperfection,  .and : 
some  possibly  unknown  orthogonal  functions  Y.  ..  The  forms  of  Y  ,  . 
Y9,  ...  are  clearly  not  harmonic  In.  x,.  but  it.  should  still  be  possible 
to  excite  mainly  the  first  characteristic:  mode  Y  using  a  lateral  1 
load  Q  applied  near  the  column  mldlength.,  The.- deflected  shape -will! 
not  be  a  half -sine  curve  as"  for  prismatic:  bars:, .  but  it  still  lwillLnat~ 
have  any  intermediate  nodes.: 

If  the  reference  level  for  laboratory  measurements  is  established 
as  before  and  the  measurements  are  taken  at.  the  point  of  load  .applica- 
tion x  =  xn,  the  measured  deflection  is  then 

QYJ(xp)  J  Y^(x„)K       1      Y3(VK1      1 

6  =  - =—  +  • + 


Kl    I  1  +  f-   Y^(x0)K2   1  +  f-    y2(x0)K3   1  ♦  |1 

X  Z.  3 


(2.16) 


E.   ANALYSIS  OF  EXPERIMENTAL  DATA 

Continuing  with  the  general  case  of  a  non-prismatic  bar,  once 
experimental  data  has  been  obtained,,  it  may  be  analyzed,  using  equation 
(2.16),  by  several  methods.   Two  specific  methods  are  outlined  here. 

Considering  the  relative  magnitudes  of  the  series  terms  in  equation 

(2.16),  the  first  is  clearly  the  dominant  term  for  any  given  value  of  T. 

If  only  the  first  term  is  retained  as  an  approximation,  6  is  a  linear 

function  of  Q  for  any  given  T  since  K  ,  Y  (x  )  and  P   are  all  constant. 

A  plot  of  6  vs.  Q,  for  a  given  T  yields  a  straight  line  with  an  inverse 

slope  equal  to 

Kl 

kp  =  l^ul  (1  +  h]  =  k  (1  +  ^}  (2-17) 

r  oJ  ri        xi 
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Thus  k      is   a  linear   function  af  Tr  indf  caring-  that:  a  plot  of:k     vs.  . T. 

p  p:- 

will  have  a  slope  equal  to  k/P  and:  an.  intercept:  equal  to  k..  By- 
simp  ly  dividing  the  intercept  by  the  slope  a  value  of  P  may  be  found. 

A  second  method,  which  may  offer  greater  accuracy,  is  a  least: 
squares  approach  using  P   as  a  parameter..  Consider  again  equation 
(2.16).   For  beams  which  are  "nearly- prismatic"  the  functions  Y  (x) 
will  be  "nearly"  sine  functions  so  that  an.  estimate  of  the  relative - 

values  of  the  first  several  coefficients  Y   (x„)K, /Y.,  (x„)K  can:. be 

n    0   1   1    On: 

made  by  assuming  the  Y  (x)  functions  are:  indeed.,  sines .   Similarly  the 

ratios  P«/P, ,  P./P, ,  pq/pT  ---»  can  be  estimated  to  be  4,.  9,.  25, . 

This  provides  a  way  of  estimating  the  contributions  of  terms  beyond  ..the 

first.   If  these  added  terms  can  be  at.  least  estimated,  they  can  be 

retained  and  used  advantageously;,  if  not,,  only  the  first  term  is? 

retained. 

For  any  given  value  of  the  parameter  P  ,  equation  (2.16)  maybe 

written 

Y2(x  )K  Y2(x  )K 

2V  (T  1     1       3V  0;  1 

+  •  =—  + 


i  +  Ji      Yi(xo)K2     L  +  !i  Y!(X0)K3  1  +  JL 

1  2  3 

Qi 
=  t-i  C.   =  Q.  p  C.  (2.18) 

k    l     11 

where  n  -  1/k.   In  the  series  C.  =  C(T.),  only  those  terms  are  retained 

which  can  be  estimated  with  some  confidence.   Let  e.  be  the  error 

i 

associated  with  the  ith  data  point. 

e.  =  Q.  p  C.  -  6-  (2.19) 

l    l  M   i    l 

Let  S  equal  the  sum  of  the  squares  of  e.  for  n  data  points. 

n 
S  =   E    e?  (2.20) 

i=l   x 
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To  find  the  most  probable   (In  the  sense  of1  least",  squares)   valuer  of fp; 
and  P      from  the  experimental  datar  proceed:  as:  fori  "lows. 

n 

T-sf  =  °  =   -z     Qt  cr  ei  (2*21) 

Z      OP  £— 1  111 

n  n 

.S,    (Q.    C.)2   P    -    .Z,    6.   Q.    C.    =  0  (2.22) 

1=1        1      1  i=l      111 

n 

£   6.    Q.    C. 
i=l   1      1     1 

P   =  "n (2:23) 

Z   (Q.    C.)2 
i=l      x      x 

Now  find  p  and  S  for  different  values  of  the  parameter  P  .   The  quantity 
S  is  a  minimum  when  the  most  probable  value  of "  F  is."  used, .  also: yield-- 
ing  a  most  probable  value  of  p.. 

For  the  case  of  a  prismatic  bar,,  equation  (2.15")  describes  the 
deflected  shape.   The  two  methods  just  outlined  are,  of  course,  applicable 
to  this  simplified  case.   All  of  the  terms  in  equation  (2.15)  are  known 
and  hence  all  may  be  used  to  improve  the  accuracy  of  the  least  squares 
analysis  if  desired. 

Both  of  the  methods  outlined  above  may  be  easily  programmed  for  the 
digital  computer.   Appendix  A  gives  program  listings  for  the  prismatic 
bar  case. 

F.   DETERMINATION  OF  THE  SECOND  BUCKLING  LOAD 

Determination  of  the  second  mode  buckling  load  using  the  theory 
outlined  above  is  straightforward  assuming  that  EI(x)  =  El(L-x).   If 
an  upward  load  -Q  is  applied  at  x  =  x   and  a  downward  load  Q  is  applied 
at  x  =  L-3C   the  resulting  equation  for  deflection  is  similar  to 
equation  (2.12).   (It  Is  convenient  to  apply  the  lateral  loads  at 
x  =  L/4  and  x  =  3L/4.)   Such  lateral  loading  does  not  excite  odd 


20 


numbered  modes,  only  evenx  numbered:  modesi,.  audi  if:  Q  is  sufficiently- 
large,  the  term  corresponding  to  the:  secondi  mode  can  be  made  : to: dominate . 

The  first,  third,  fifth,.  terms  are  negligible,  and  the  fourth,  .sixth, 

...  are  small  compared  to  the  second..  Thus,,  the  sort  of  data  reduction 
discussed  above  will  lead  to  an  evaluation  of:  the  second  mode-.. 

If  EI(x)  i-   EI(L-x)T  that  is  if  there  is:  no:  symmetry  with  respect  ~to 
a  normal  plane  at  x  =  L/2r  it  may- he,  quite:  difficult  to  apply  adateral  _ 
loading  system  which  does  not  excite  the  first:  mode .   A  suggested :ap^ - 
proach  is  to  apply  an  upward  load  at  x  —  L7M--  and  an  equal  downward  :load 
at  a  point  "near"  x  =  3L/M-..   Unless  the  second,  point  is  selected: 
correctly,  the  mixture  of  first:  and  secondf  modes  thus  excited; will  '. 
not  provide  a  linear  plot  of  data..  One  would,"  have  to  try  several  .such 
points  in  order  to  obtain  the  desired  sort  off  plot  which  could  then 
be  interpreted  so  as  to  yield  the  second  buckling  load. 
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III.      EXPERIMENTAL  PEOCEDURF: 

A.  DESIGN 

During  the  design  and  fabrication,  of  the  test:  apparatus ,  strict 
attention  was  paid  to  simplicity..  The  major  problems:  which  arose 
where  how  to  best  simulate  pinned  end  conditions:  ihr.  the  laboratory,  . 
and  how  to  apply  the  required  axial  and  laterall  loads:. .  Several-designs: 
were  considered  which  could  easily  transmit  an  applied:,  compressive. 
axial  load,  but  they  could  not  easily-  support  an  applied  lateral  load: 
without  introducing  unknown  end  moments..  Application  of  a  compressive 
axial  load  using  a  universal  testing  machine,  could:  be:  accomplished,  . 
but  lateral  loads  could  not  be  easily  applied  or  measured.   Dead  weights 
seemed  to  be  the  simplest  means  of  applying  both  axial  and  lateral 
loads.   The  final  design  of  the  test  apparatus  Is  shown  schematically 
in  Fig.  (3.1a). 

B.  FABRICATION 

Compressive  axial  loads  could  not  be  easily  applied  using  dead 
weights  so  it  was  decided  to  use  tensile  loading  for  convenience.   This 
also  allowed  determination  of  higher  mode  buckling  loads.   However, 
there  is  a  penalty  for  using  tensile  axial  loads.   This  will  be  discussed 
in  Section  IV- B.  Pinned  end  conditions  were  simulated  using  steel 
cylinders  on  either  end.   The  cylinders  were  polished  and  rested  on 
a  finished  surface.   A  mixture  of  graphite,  molybdenum  sulfide,  and 
machine  oil  was  used  as  a  lubricant  to  eliminate  dry  friction  between 
these  cylinders  and  the  test  bed..   Flexible  steel  cable  was  used  to 
transmit  the  axial  load  T.   The  cable  was  brazed  to  a  steel  ring  for 
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attachment  to  the  end  cylinders .   A.  slat  was  cut  itr  each-  exrdi  cylinder., . 
perpendicular  to  the  longitudinal  axis  and  the  ring  was  inserted:..  A 
pin  which  provided  a  "knife  edge"'  contact  was  inserted  parallel!  to:  the-: 
longitudinal  axis  of  the  cylinder,  passing  through  the  eye  off  the 
ring,  see  Fig.  (3. lb).   The  knife  edge  of  the  pin  ran  directly  along 
the  longitudinal  axis  of  the  cylinder.   The  effective  length-  of  the: 
column  was  measured  from  centerline  to  centerline  off  the:  end  cylinders 
and  was  19.00  inches.   Since  a  near-zero  moment  exists  near:  these: 
centers,  changes  in  EI  near  the  ends  have  negligible  effect:  on  he ndihgj . 
and  the  effective  column  length  may  be  safely  considered  as  the  distance 
between  centers  of  the  end  cylinders.   The  column,  itself  was  1/4-""  x.  1/8" 
cold  drawn  1018  steel,  with  a  modulus  of  elasticity  E~  —  23". .4-  x  ID6'  psi.. 
It  was  brazed  to  the  end  cylinders.   The  lateral  load  was  applied  with 
a  knife  edge  and  dead  weights  as  shown  in  the  figure.   Deflection 
measurements  were  taken  with  a  simple  wire  probe  attached  to  a  verti- 
cally mobile  micrometer.   The  knife  edge  assembly  used  for  application 
of  lateral  loads  had  a  hole  drilled  in  the  flat  upper  surface  which  was 
filled  with  mercury.   The  probe  was  lowered  until  contact  with  the 
mercury  was  made,  completing  a  simple  neon  light  bulb  circuit,,  see 
Fig.  (3.2).   For  second  mode  buckling  tests,  two  lateral  loads  were 
applied,  one  at  each  quarter  point.   The  downward  load  was  applied  as 
before.   The  upward  load  was  also  applied  using  a  knife  edge  and  dead 
weights.   A  pulley  was  used  to  convert  dead  weight  into  upward  force. 
Since  an  upward  force  was  applied  at  one  end  of  the  beam,  the  cylinder 
on  that  end  had  to  be  restrained  against  upward  motion.   This  was  done 
with  a  simple  "hood",  see  Fig.  (3.3). 
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Figure  3.2 
APPARATUS  FOR  DETERMINATION  OF  FIRST   BTJ3KLING  LOAD 
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Figure  3.3 
APPARATUS  FOR  DETERMINATION  OF  SECOND  BUCKLING  LOAD 
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C.   TESTS 

The  reference  level  for  testing  was  established  before  each  test 
by  first  applying  the  axial  load  T,  and  then  taking"  a  reading  from 
the  micrometer  with  the  knife  edge  assembly  in.  place.   The  weight 
of  the  knife  edge  caused  some  small  deflection,,  but  this  may  be 
treated  in  the  same  manner  as  the  initial  imperfections  were.   Once 
the  reference  level  was  established,  the  lateral  load  was  increased 
from  0.0  to  4.0  pounds  in  half-pound  increments  with  a  final  1/4 
pound  increment  to  give  a  maximum  load  of  4-25  pounds.   No  greater 
lateral  load  was  applied  so  as  to  assure  that  plastic  bending  did 
not  take  place.   This  was  done  for  15  different  values  of  T  ranging 
from  0.0  to  42.0  pounds,  in  three-pound  increments.   For  each  such 
loading,  the  values  of  loads  and  corresponding  micrometer  readings 
were  recorded. 

The  same  method  of  establishing  a  reference  was  used  for  second 
mode  buckling  tests.   An  upward  load  (-Q)  was  applied  at  x  =  L/4  and 
a  downward  load  (Q)  was  applied  at  x  =  3L/4.   The  value  of  Q  ranged 
from  0.0  to  6.5  pounds,  in  half-pound  increments.   The  range  of  Q 
could  have  been  extended  another  two  pounds  without  causing  plasticity, 
but  the  effects  of  friction  became  evident  above  6.5  pounds,  since 
increases  in  Q  caused  almost  immeasurable  changes  in  deflection.   The 
values  of  T  ranged  from  0.0  to  41.0  pounds  in  five-pound  increments, 
the  odd  pound  arising  from  the  weight  pan  used.   A  few  anomalies  in 
the  data  for  second  mode  were  noted  and  will  be  discussed  subsequently. 
Several  attempts  to  find  the  cause  of  the  discrepancies  were  made,  but 
proved  fruitless. 
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Once  the  data  was  ohtained:,.  it:  was:  analyzed: using  both  of  r  the-: 
methods  described  previously.. 

As  a  means  of  substantiating;  earlier  test:  results ,  another,  type 
of  experiment  was  done  to  determine,  the  primary  buckling  load,  using 
axial  comDression  and  interpreting  the  data:  according  to  Southwell'.s 
method.   The  apparatus  used  for  this  test:  is:  shown;  in  Fig.  (3.4).   A 
cantilever  beam,  embedded  in  a  vertically  mobile:  "wall"  was  used  to  : 
apply  a  compressive  axial  load  to  the'  column..  The  cantilever: was: 
instrumented  with  four  strain-  gages.,,  two:  mounted"  parallel  to  the 
direction  of  principal  strain,,  and  two  perpendicular  to  it  to  compen-  - 
sate  for  Poisson''s  effect..  The  gages  were:  connected  in  a  four-arm  - 
bridge  circuit  and  hooked  to  a  Baldwin  Model"..  SR-4-- Strain  Indicator: -- 
Type  N.   The  column  was  also  instrumented  with  strain  gages.   One  gage 
was  mounted  on  each  side  of  the  column  at  mid-span,  parallel  to  the 
direction  of  principal  strain.   These  two  gages  were  designed  specif- 
ically for  use  on  mild  steels  and  were  self-compensating  for  small 
ambient  temperature  changes.   These  two  gages  were  connected  in  a  two- 
external-arm  bridge  hooked  to  an  Ellis  Bridge  Amplifier  and  Meter 
(BAM-1).   The  output  of  the  BAM  was  monitored  on  a  digital  voltmeter. 

The  cantilever  strain  gages  were  calibrated  using  dead  weight  tip 
loads,  providing  a  direct  correlation  between  indicated  strain  and 
applied  load. 

Before  each  test,  the  no-load  reference  for  the  cantilever  bridge 
was  noted.   A  small  load  was  then  applied  to  hold  the  column  vertical, 
and  the  reference  reading  for  the  column  bridge  was  noted.   Preliminary 
tests  indicated  that  the  Southwell  plot  became  linear  for  P  >_  .  5P  . 
Therefore,  the  axial  load,  including  the  weight  of  the  end  cylinder, 
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was  increased  from  16.  to  26.  pounds  in  one-pound  increments.   At 
each  load  increment  the  strain  readings  from  the  cantilever  bridge 
and  the  voltage  readings  from  the  column  bridge  were  recorded,  the 
strain  reading  being  proportional  to  applied  axial  load,  and  the 
voltage  change  (AV  =  V  -  V  )  being  proportional  to  column  deflection. 

AV/P  versus  AV  was  then  plotted  for  each  of  four  runs  and  the 
value  of  P  was  determined  by  the  reciprocal  of  the  slope.   A  computer 
program  was  utilized  to  obtain  the  best  fitted  line,  for  each  run, 
again  using  a  linear  least  squares  fit,  and  the  average  value  of  P 
for  all  runs  was  calculated  (see  Appendix  A). 

A  standard  tensile  test  was  run  on  a  specimen  of  the  material 
used  for  the  column  to  determine  Young's  modulus  experimentally.  The 
result  was  E  =  (29.4  +  .1)  x  106  psi. 
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IV.   RESULTS  ANT  DISCUSS  ION" 

A.   RESULTS 

The  figures  tabulated  In  Table  I  represent  the  experimentally 
determined  values  of  the  buckling  load  for  both  first  and  seconi 
modes.   As  described  earlier  (section  2-E)  two  methods  were  used.ro: 
analyze  the  data:   one  being  a  linear  approximation,,  frereafter 
method  I;  the  other  being  a  least  squares  analysis  rerainihg:  any- 
desired  number  of  terms  In  equation  (2.15),.  hereafter  method  III.. 
The  average  value  of  P  arrived  at  by  use  of  Southwell's  method: is 
also  tabulated. 


Southwell ' s 
Method  I      Method  II         Method 

Number  of 

terms  retained  1  14  


1st  mode 

buckling  load  (lbs.)     31.0       32.0     31.2       32.1* 

2nd  mode 

buckling  load  (lbs.)    127.8      110.5    109.0        


Table  I 
Summary  of  Results 

Typical  plots  of  6  vs.  Q  and  k  vs .  T  for  method  I  are  shown  in  Fig. 

p  & 

(4.1)  and  (4.2)  respectively.  A  typical  plot  for  determination  of  P 
by  Southwell's  method  is  shown  in  Fig.  (4.3).  Compiled  data  for  all 
tests  is  listed  in  TablesII  -  V. 


Average  value  of  four  tests, 
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E&.TS.  FOR  SOUTHWELL' 'ST  METHOII 


CSR 


V 


AV   AV/P 


CSR 


F 


V 


AV 


AV/P 


1226 

16.3 

-  55 

49 

3..0Q 

1226 

16.3. 

-  61 

53 

3.25 

1211 

17.3 

-  64 

57 

3.29 

1226 

16.3 

-  66 

58: 

3.35 

1197 

18.3 

-  72 

65 

3. .55 

1197 

18.3 

-  73 

65 

3.55 

1182 

19.3 

-  84 

77 

3..99 

1182 

19. 3 . 

-  81 

73 

3.78 

1168 

20.3 

-  94 

87 

4.29 

1168 

20.3 

-  92 

84 

4*14 

1153 

21.3 

-1Q5 

98 

4- ..60 

1153 

21.3 

-106 

98 

4,60 

1139 

22.3 

-126 

119 

5. .33 

1X3291 

22.3 

-123 

115 

5.16 

1124 

23.3 

-144 

137 

5..8ES 

ri2H- 

23.3 

-140 

132: 

5.67 

1110 

24.3 

-163 

156 

6.42 

1110 

24.3: 

-162 

154 

6.34 

1095 

25.3 

-186 

179 

7..08 

1095 

25,  3": 

-197 

191. 

7755 

1081 

26.3 

V   = 
o 

-224 
-7mv 

217 

8..  25 

1081 

25.3- 

V  = 
o: 

-227 
=  -8mv 

219 

8.33 

CSR 


V 


AV 


AV/P 


CSR 


P 


V 


AV 


AV/P 


1226 

16.3 

-  46 

45 

Z..76 

1226 

16.3 

-  55 

50 

3.07 

1211 

17.3 

-  54 

55 

3.12 

1211 

17.3 

-  61 

56 

3.24 

1197 

18.3 

-  65 

64 

3.50 

1197 

18.3 

-  69 

64 

3.50 

1182 

19.3 

-  74 

73 

3.78 

1182 

19.3 

-  80 

75 

3.89 

1168 

20.3 

-  96 

95 

4.68 

1168 

20.3 

-  89 

84 

4.14 

1153 

21.3 

-106 

105 

4.93 

115  3 

21.3 

-106 

101 

4.74 

1139 

22.3 

-116 

115 

5.16 

113  8 

22.3 

-121 

116 

5.20 

1124 

23.3 

-143 

142 

6.09 

1124 

23.3 

-138 

133 

5.71 

1110 

24.3 

-159 

158 

6.50 

1110 

24.3 

-160 

155 

6.38 

1095 

25.3 

-180 

179 

7.08 

1095 

25.3 

-181 

176 

6.96 

1081 

26.3 

V   = 
o 

-209 
-Imv 

208 

7.91 

1081 

26.3 

V   = 
o 

-211 
-5mv 

206 

7.83 

Note:   CSR  =  Cantilever  strain  reading  (yin/in) 
P  =  Compressive  axial  load  (lbs) 
V  =  Voltage  reading  from  digital  voltmeter  (mv) 

An  empirical  equation  based  on  earlier  calibration  tests  was  used  to 
relate  CSR  and  P 

P(lbs)  =  (CSR  -CSR)xl.01b/14.76uin/in+Wt 
o  c 

For  all  tests  CSR  =  1440uin/in. 

Wt   =  weight  of  the  end  cylinder  =1.5  lbs. 


AV  =    V-V 


Table  IV 
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DATA  FOR  TENSILE  TEST  TQ 
DETERMINE  YOUNG'S  MODDLOS 

SRIL     SRDL      T       SIL      SDL      ASPd    STRESS. 


0 

1920 

1922 

0 

0 

0 

0 

o: 

100 

1700 

1700 

100 

220 

222 

m.s 

3226 

200 

1478 

1482 

200 

442 

440 

Z2D...5 

6452. 

300 

1255 

1262 

300 

665 

660 

33TE..3: 

9677" 

400 

1038 

1043 

4Q0 

882 

S7"9 

440.  .3": 

12903: 

500 

819 

827 

500 

HOI 

IG35 

549..  0' 

1B129 

600 

602 

610 

600 

1318 

1312 

65.T.5 

19355i 

700 

382 

385 

700 

1538 

1537 

76  &.  & 

225  BE 

800 

161 

163 

8QQ 

1759 

1759 

879.. 5 

25806 

850 

53 

53 

850 

1867 

1869 

934,. 0: 

27419 

Note:   Strain  readings  (uin/ln)  were  taken  with  increasing;  EoacL 
(SRIL)  and  decreasing  load  (SRDL).. 

T  =  Tensile  axial  load  (lbs) 
ASPG  =  Average  strain  per  gage 

The  cross-sectional  area  of  the  specimen  was  0.031  in  . 

SIL  =  Actual  strain,  increasing  load 

SDL  =  Actual  strain,  decreasing  load 


Table  V 
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B.   DISCUSSION 

In  order  to.  establish  a.  point  of"  CDraparikorr.  between  experimental 
and  theoretical  results,,  the  Euler-  load,  may  he:  calculated,  based  on: 
the   following  measured  quantities : 

E   =  Young's   Modulus   =  ( 29  ..4  +^  ..I)    x.  10€   psi 

w  =  column  width   =  0..2.4SF  -h  ..OOTl  in.. 

h  =  column  depth   -   Q..1Z5   i-_  ..001  in.. 

L  =  effective  column.  length.  =  I3..00'.  -f-  .Dll  in.. 


pi  = 


tt^EI        ttzE      wh 


2      12 


=  3Z..44-  lbs-. 


(4.1) 


L^  L' 

Due   to  the  uncertainties    In.  the  measured:  quantities  ,   there   is  ;  an: 
uncertainty  In  the   calculated  value-  of  F  ...     The  magnitude   of.  this, 
uncertainty  may  be   approximated  using  a  "logarithmic"   error  analysis 
Taking  the   logarithmic  derivative   of  equation    (4.1) 


dpll 

dE 
E 

+ 

dw 
w 

+   3 

dh 
h 

+      2 

dL 
L 

(4.2) 


dP      =   P    (.0324)    =    1.05    lbs 
P      =    32.44   +   1.05    lbs. 


(4.3) 


The  experimental  values  of  P   as  determined  by  methods  I  and  II  fall 
only  slightly  outside  the  range  of  uncertainty  associated  with  the 
calculated  Euler  load.   The  value  of  P  determined  by  Southwell's  method 
may  be  considered  as  an  upper  bound  in  view  of  Southwell's  theory  [5] 
which  indicates  experimental  values  will  be  either  exact  or  high. 

The  Euler  load  for  second  mode  buckling  is  four  times  that  for 
first  mode.   Thus,  we  would  expect 

P   =  129.32  +_  4.20  lbs. 
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The  experimental  value  determined:  by  method  TJ  is:  well-  within:  the: 
range  of  uncertainty*  hut.  that:  determined,  hy  method:  II  seems  "con- 
siderably In  error.   This  discrepancy  is:  attriiiutable  to  inconsis-- 
tencles  in  measured  deflection..  Referring-  to  Table  III,  particularly 
the  data  for  T  =  11..Q  and  T  =  1E..0  lbs..,,  it  may  be  seen  that  increas- 
ing T  resulted  In  Increased-  deflection..  For:  every  given  value  of:' 
lateral  load  Q,  the  deflections-  measured  with  T  =  16.0  lbs.  are 
approximately  -G03  Inches  greater  than  those:  measured  for  T  =r111.0' 
lbs.   A  second  test  with.  T  —  16".  .0'  Us.,  was:  run  on  another  day  as:a: 
check,  and  the  values  obtained  earlier  were-  repeated.   This  discrepancy 
was  completely  unexpected  In  view  of/ the  theory,,  so  several  checks  of. 
the  experimental  apparatus  were  made,,  but  no:  unexpected  sources  of: 
error  came  to  light..  This  inconsistency  occurred  for  several- (but 
not  all)  values  of  increased  T. 

Despite  these  anomalies,  the  data  can  be  treated  using  method  I. 
For  a  given  value  of  axial  load  (T),  a  plot  of  deflection  (6)  vs. 
lateral  load  (Q)  is  made  but  only  the  slope  of  this  line  is  considered. 
For  this  reason  a  constant  error  in  measured  deflection  has  no  effect 
on  the  outcome  of  the  analysis..  Method  II,.  on  the  other  hand,  relies 
on  the  absolute  rather  than  the  relative  magnitudes  of  measured 
deflections,  and  gives  erroneous  values  for  P  ,  as  expected. 

The  theory  developed  earlier,  particularly  equations  (2.13)  and 
(2.1M-),  shows  why  it  is  advisable  to  retain  several  series  terms  when 
the  axial  loading  is  tensile.   For  any  given  tensile  loading,  the 
effect  of  additional  series  terms  in  equation  (2.14)  is  greater  than 
that  for  compressive  loading  (Eq.  2.13).   Thus  a  single  term  approxi- 
mation may  be  more  easily  justified  for  compressive  loads  than  for 
tensile. 

42 


V.   CONCLUSIONS  AND  RCCDMMQIITATTONS' 

Experimental  determination  off  the  primary  buckling'  load-off  a-:bar>' 
can  be  accomplished  without  reliance  on  some  accidentai,  initial 
imperfection.   Lateral  loading  is  used  to:  assure  dominance  of:  the - 
primary  mode  in  the  deflected  shape..  The  axial  load  may  be  either: 
tensile  or  compressive-   By  correlating:  observed,  deflection-  and : 
applied  lateral  load  for  several  distinct"  values  off  axiall  load,  .  the: 
buckling  load  may  be  easily  and  accurately  determined.. 

If  a  tensile  axial  load  is  used,,  experimental,  determination  of  : 
the  second  buckling  load  is  possible,,  using"  appropriate  lateral,  loads 
to  assure  dominance  of  the  second  mode  in  the  deflected  shape. 

The  developed  theory  and  experimental  technique,  are  generally 
applicable  to  bars  of  any  cross-section,,  although  orthogonal  load- 
ing schemes  for  determination  of  higher  mode  buckling  loads  may  be 
difficult  experimentally  for  non-prismatic  bars.   The  use  of  compress- 
ive  axial  loading  is  advantageous  for  non-prismatic  bars  since  it 
minimizes  the  effects  of  possibly  unknown  higher  modes  in  the  deflected 
shape. 

It  appears  that  techniques  suggested  may  also  be  extended  for  use 
in  determining  plate  buckling  loads. 
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APPENDIX  A  -  COMPUTES  PROGRAM'  LISTINGS: 

The  computer  programs  listed  on  the  following  pages  were  used. 
to  analyze  experimental  data..  All  programs,  are  written  in  FORTRAN 
IV  for  use  on  an  IBM  3 60  computer. 

The  programs  for  Method  I  and  Southwell's:  Method  are  essentially 
the  same.   The  program  was  taken  from  McGalTa.  [33_, .  in  part,,  with  ae 
few  additions  and  minor  changes  made  as  necessary.. 

The  input  required  for  each,  program,,  excluding  data,  is  minimal - 
and  is  explained  In  the  body  of  each  program..  (Note  that  in  the: 
programs  the  symbol  P  Is  used  to  denote  axial- tension. ) 
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APPENDIX  B  -  EVALUATION:  QF  INTEGRALS  INVOLVING  ORTHOGONAL^ 

FUNCTIONS  AND  THEIR  DERIVATIVES: 

The  equation  and  conditions  governing-  the.  shape:  off  the--  strained  : 
centerline  of  a  pin-ended  column,  subjected  ta  m  axial  compressive: 
load  P  are 

BY"  +  PY  =  Q  where  B  =  B(x>  =  EI(x)  (Br  la  I 

Y(0)  =  Y(L)  =  0  (Brlb) 

Let   (Y.  ,   P.)   and  (Y.,   P..)   be  two   distinct   solutions   to   this:  equation- 
113        3 

BY."'   +   P..Y..    =  0 
1  11 

BY."    +   P.Y.    =   0  (Br2) 

Multiplying  both  sides  of  equation  B-2  by  Y." 
L  L 

/  BY."  Y."  dx  +  /  P.Y.Y."  dx  =  0 
J  1    3       J       1  1  3 

0  0 

L  L  L 

P.  /   Y.Y."  dx  =  f   BY."Y."  dx  =  -P.  f    Y."Y.  dx         (B-3) 
J    1  3         3j    1   3 

0  0  0 


-L  -L 


.  /   Y.Y." 
ij     1  3 

.  f  Y."Y. 
1  I         13 

Fy.y. 
L  3  1 


.'.  P.  /   Y.Y."  dx  =  P_.  /   Y_."Y,  dx  (B-4) 

J. 

=  P.,' 


-  /   Y . ' Y . '  dx 

0    J     x   ' 


=  -P.  f     Y.'Y.*  dx 

"I     *   : 

Also  the  left  side  of  equation  B-4  is 

L  L 

P.  /   Y.Y."  dx  =  -P.  f     Y.'Y.*  dx 

.'.  (P.  -  P.)/   Y.'Y.'  dx  =  0 
1     3J    1   3 

0 
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But:  P_  and:  F..  atre  dxs1ri±rc±-_  (ii  i=  f) 


..    /    v--  'Y'  '' 
I        i'     I 

ATsor   from  equation  B-4- 
L  L. 


_  ,     Y.,' V  dx   -  o:       ii  £  j 


f   Y."Y..  dx   =-/       Y."Y..  dx :=  0  i   ?M 

J       x      j  J  x.     j 


Q 
By  equal; ion.  H— 3 
L 


0 


[    BY.MY."    dx:  =  -P../ 


dx:  =  -P.../        Y.."Y...  dx  =   0        i. /=  j 
I  ^    I 


Q  o: 

Summarizing  the  proceedings  arguments 


(B-5) 


(B-6) 


(B-7) 


i  Y.Y."  dx  =  I. .6.....;:  I\.  =/  Y.Y."  dx 

/  Y.'Y.1  dx  =  J.  6.  .;  J.  =■/  (Y.')2dx 

J  i  D        x  xj   x  J  i 

0  0 


BY."Y."dx  =  K.6..;  K.  =/   (Y.")2dx 
1   ]       l  13    l  1 


where   6 . .   is  the  Kronecker  delta 
1: 


13 


V 


0  i  t   j 

1  x  =  j 


Returning  to  equation  B-3  and  B-4  it  may  be  seen  that 


K.  =  P.J. 
1     11 
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